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ABSTRACT 

We present time-resolved spectroscopy and photometry of the double-lined eclipsing 
cataclysmic variable V347 Pup (= LB 1800). There is evidence of irradiation on the 
inner hemisphere of the secondary star, which we correct for using a model to give 
a secondary star radial velocity of Kr — 198 ± 5 kms^^. The rotational velocity 
of the secondary star in V347 Pup is found to be usini = 131 ± 5 kms~^ and the 
system inclination is i = 84.0° ± 2.3°. From these parameters we obtain masses of 
Ml = 0.63 ± 0.04 Mq for the white dwarf primary and Ma = 0.52 ± 0.06 Mq for 
the M0.5V secondary star, giving a mass ratio oi q — 0.83 ± 0.05. On the basis of 
the component masses, and the spectral type and radius of the secondary star in 
V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows 
many of the characteristics of the SW Sex stars, exhibiting single-peaked emission 
lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. 
We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to 
be close to the maximum allowed in a pressureless disc. 

Key words: accretion, accretion discs - binaries: eclipsing - binaries: spectroscopic 
- stars: individual: V347 Pup - novae, cataclysmic variables. 



1 INTRODUCTION 

Cataclysmic variables (CVs) are close binary stars consist- 
ing of a red-dwarf secondary transferring material onto a 
white dwarf primary via an accretion disc or magnetic ac- 
cretion stream. V347 Pup is an example of a nova-like vari- 
able (NL), a class of CV with high mass transfer rates 
and no recorded nova or dwarf-nova type outbursts; see 
[Warner (1995a)|for a comprehensive review of CVs. 
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A knowledge of the masses of the component stars in 
CVs is fundamental to our understanding of the origin, 
evolution and behaviour of these systems. Population 
synthesis models (e.g. |Kolb, King fc Baraffe 2001| | and the 
disrupted magnetic braking model of CV evolution (e.g. 
[Spruit fc Hitter 1983| Rappa port, Verbunt fc Joss 1983| 
can be observationally tested only if the number of reliably 
known CV masses increases. One of the most reliable ways 
to measure the masses of CVs is to use the radial velocity 
and the rotational broadening of the secondary star in 
eclipsing systems. The radial velocity of the disc emission 
lines is often an unreliable indicator of the white dwarf 



© 0000 RAS 



2 T. D. Thoroughgood et al. 



motion because of contamination from, for example, the 
bright spot. At present, reliable masses are known for only 
^20 CVs, partly due to the difBculties in measurement (see 
IHim th & Dhi llon 19981 for a review). 

V347 Pup was identified spectroscopically as a NL by 
[Buckley et al. (1990)| from the presence of high-excitation 
emission lines. Even though V347 Pup emits at X-ray wave- 
lengths (as the Uhuru X-ray source 4U 0608-49), the NL 
classification was favoured over a magnetic CV class on ac- 
count of the negligible polarisation present. The study by 
[Buckley et al. (1990)| revealed a bright and deeply eclips- 
ing system, with a spectroscopic and photometric orbital 
period of 5.57 hrs. Their measured system inclination and 
emission-line radial velocity curve, together with an empir- 
ical secondary star mass estimated from the orbital period, 
suggested a high primary mass close to the Chandrasekhar 
limit. 

A multiwavelength study by [Mauche et al. (1994)[ re- 
vealed an X-ray spectral energy distribution similar to many 
dwarf novae in outburst, with a likely origin in an ex- 
tended emission region rather than the boundary layer. The 
UV emission lines appear to have a similar origin and, in 
a later paper by [Shlosman, Vitello fc Mauche (1996)[ their 
behaviour in eclipse was sucessfuUy modelled as disc light 
scattered in a rotating wind. The presence of an accre- 
tion disc in V347 Pup was confirmed by a rotational dis- 
turbance of the optical emission lines through primary 
eclipse UMauche et al. 19941 [Stih, Buckley fc Garlick T998l l. 
The latter authors found evidence for spiral arms and disc 
overflow accretion, and identified the low-excitation optical 
emission profiles as a composite of emission from the accre- 
tion disc and secondary star. 

Secondary star absorption lines were found by 
[Diaz fc Hubeny (1999)[ who measured the system pa- 
rameters of V347 Pup using the radial velocity semi- 
amplitudes of the primary and secondary stars. The ra- 
dial velocity of the optical emission lines in V347 Pup 
varies widely in the literature, with published values 
of 134 ± 9 kms"^ fB uckley et al. 1990| l, 122 ± 19 
kms"^ Pvlauchc ct al. 1994}, 156 ± 10 kms'^, 125 ± 13 
kms~^ fStiU, Buckley fc Garlick 1998t and 193 ± 16 kms"^ 
l |Diaz fc Hubeny 1999|l. The radial velocities of the UV emis- 
sion lines published by [Mauche et al. (1994) [ ranged between 
220 and 370 kms~^ with large phase shifts between spectro- 
scopic conjuction and photometric mid-eclipse. This wide 
range in values, and the known unreliability of using disc 
emission lines in NLs to determine the motion of the white 
dwarf (e.g. iPhill on, Marsh fc Jones 1997| l, makes the deter- 
mination of system parameters from the secondary star fea- 
tures alone highly desirable. In this paper, we derive the 
system parameters from the radial and rotational velocities 
of the secondary star in V347 Pup. 



2 OBSERVATIONS AND REDUCTION 

During January and December 1998 and January 1999, we 
obtained optical spectra of V347 Pup using the Cassegrain 
spectrograph -|- SITel CCD chip on the SAAO 1.9-m tele- 
scope. Simultaneous photometry was available for most of 
the spectra using the SAAO 1.0-m telescope with the TEK8 
CCD chip. See Tableland its caption for full details. 



On the December 1998 run, we observed 17 spectral 
type templates ranging from G7V-M5.5V and telluric stars 
to remove atmospheric features. We observed flux standards 
on both the 1.9-m and 1.0-m telescopes on all nights. 

The spectra and images were reduced using stan- 
dard procedures (e.g. [Dhillon, Jones fc Marsh 1994| 
[Thoroughgood et al. 200T] |. The photometry data were 
corrected for the effects of atmospheric extinction by 
subtracting the magnitude of a nearby comparison star. 
The absolute photometry is accurate to approximately ±0.5 
mjy; the relative photometry to ±0.01 mag. Comparison 
arc spectra were taken every ^^40 min in order to calibrate 
the wavelength scale and instrumental flexure. The arcs 
were fltted with fourth-order polynomials with an rms 
scatter of better than 0.04A. Where possible, slit losses were 
then corrected for by multiplying each V347 Pup spectrum 
by the ratio of the flux in the spectrum (over the whole 
spectral range) to the corresponding photometric flux. 



3 RESULTS 
3.1 Ephemeris 

The times of mid-eclipse for V347 Pup were determined by 
fitting a parabola to the eclipse minima in the photometry 
data. A least-squares fit to the 21 eclipse timings listed in 
Table |5| yields the ephemeris: 



Tmid-ociipso = HJD 2 446 836.96176 ± 0.231936060 E 
± 0.00009 ± 0.000000006. 



(1) 



Our new ephemeris is exactly the same as that given by 
[Baptista fc Cieslinski (1991)[ except we have reduced the 
errors on both the zero point and orbital period. We find 
no evidence for any systematic variation in the O-C values 
listed in Tabled 

3.2 Average spectrum 

The average spectra of V347 Pup, uncorrected for orbital 
motion, are shown in Fig.0 In Table[31 we list fluxes, equiv- 
alent widths and velocity widths of the most prominent lines 
measured from the average spectra. 

The Balmer emission lines are broad, symmetric and 
single-peaked, instead of the double-peaked proflle one 
would expect from a high inclination accretion disc (e.g. 
IHorne fc Marsh 198611 . This behaviour is characteristic 
of the SW Sex stars (e.g. [Dhillon, Marsh fc Jones 1997ll. 
Previ ous studies of V347 P up by [Buckley et al. (1990)] 
and [Diaz fc Hubeny (1999)[ agree with this single- 
peaked observation, however, the study by 
[Still, Buckley fc Garlick (1998)[ shows double-peaked 
low excitation lines (although this could be due to the 
presence of absorption cores) . The Hel A6678A line appears 
to be composed of a narrow single-peaked component 
superimposed upon a broad double-peaked component. 
The other Hel emission lines can clearly be seen in the 
wavelength region centred on A4610A as double peaked in 
nature, with the possible exception of Hel A447lA. High- 
excitation emission is present through Hel A4686A and the 
CIII/NIII AA4640-4650A Bowen fluoresence complex. 

The secondary star is clearly visible in the aver- 
age spectra as absorption lines of the neutral metals 
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Table 1. Journal of observations. During Jan 1998, we used Grating No. 4 to give a wavelength range of ~4200-5060A (Xcen = 4610A) 
at 0.99-A (64 kms~^) resolution. Grating No. 4 was again used on 28 Dec 1998 and 23 Jan 1999 to give a wavelength range of ^4900— 
5720A (Xcen = 5290A) at 0.95-A (54 kms"!) resolution. On 25 and 27 Dec 1998, we used Grating No. 5 to give a wavelength range of 
~5960-6725A (Acen = 6330A) at 0.88-A (42 kms~^) resolution. Simultaneous photometry for the December 1998 spectra was recorded 
in the Johnson-Cousins V and R bands. Photometry was also available during the January 1998 run in the Stromgren b and y filters. 
The seeing measured around 1.0 arcsec with photometric conditions on 25 and 27 Dec 1998 and Jan 23 1999. On 28 Dec 1998, however, 
the seeing was poor and patchy high cloud was present. The seeing varied between 1.0—1.5 arcsec over the January 1998 run. The epochs 
are calculated using the new ephemeris presented in this paper f equation lit. 



UT Date 1.9-m No. of Exposure Epoch Epoch 1.0-m No. of Exposure Epoch Epoch 

Acen (A) spectra time (s) start end filter images time (s) start end 



1998 Jan 07 


4610 


113 


200 


17178.62 


17179.91 


b 


843 


30 


17178.88 


17179.93 


1998 Jan 08 


4610 


32 


200 


17182.88 


17183.23 


b 


105 


30 


17183.83 


17184.04 


1998 Jan 10 


4610 


117 


200 


17191.94 


17192.94 


y 


1379 


30 


17191.61 


17192.96 


1998 Jan 11 


4610 


129 


200 


17195.74 


17197.25 


b 


1217 


30 


17195.81 


17197.08 


1998 Jan 12 


4610 


116 


200 


17200.09 


17201.44 


b 


373 


30 


17200.14 


17200.52 


1998 Dec 25 


6330 


68 


300 


18696.47 


18697.55 


R 


962 


30 


18696.32 


18697.65 


1998 Dec 27 


6330 


61 


300 


18705.18 


18706.17 


R 


748 


30 


18704.95 


18706.25 


1998 Dec 28 


5290 


23 


300 


18709.63 


18710.16 


V 


442 


30 


18709.50 


18710.29 


1999 Jan 23 


5290 


64 


300 


18821.30 


18822.49 




no 


photometry available 





Table 2. Times of mid-eclipse for V347 Pup according to Buckley 
et al. (1990; B90), Baptista & CiesUnski (1991; BC91) and this 



paper. 


Cycle 


HJD 


Uncertainty 


O-C 


Reference 


(E) 


at mid— eclipse 
(2,400,000+) 


on HJD 


(sees) 




-4 


46836.0379 


5x10-"* 


335.96 


B90 





46836.9621 


5x10--* 


29.74 


B90 


39 


46846.0059 


5x10""' 


-117.69 


B90 


43 


46846.9333 


5x10-** 


-147.43 


B90 


48 


46848.0930 


5x10^"* 


-145.73 


B90 


56 


46849.9500 


5x10--* 


-15.13 


B90 


65 


46852.0373 


5x10--* 


-25.89 


B90 


69 


46852.9651 


5x10-* 


-21.08 


B90 


78 


46855.0533 


5x10-* 


45.92 


B90 


6177 


48269.63136 


1.5x10--* 


48.28 


BC91 


7583 


48595.73325 


1.1x10-* 


30.04 


BC91 


7587 


48596.66022 


9x10-5 


-36.85 


BC91 


17179 


50821.39199 


5x10-* 


56.29 


This paper 


17184 


50822.55127 


5x10-* 


21.70 


This paper 


17192 


50824.40676 


5x10-* 


21.83 


This paper 


17196 


50825.33428 


5x10-* 


2.46 


This paper 


17197 


50825.56611 


5x10-* 


-6.70 


This paper 


18697 


51173.47035 


1x10-* 


6.20 


This paper 


18705 


51175.32562 


1x10-* 


-12.68 


This paper 


18706 


51175.55802 


1x10-* 


26.54 


This paper 


18710 


51176.48519 


1x10-* 


-22.21 


This paper 



Cal, Fel and Mgl, as seen in [Diaz fc Hubeny (1999)| Sec- 
ondary star features in the SW Sex stars are not un- 
usual in the longer period systems, such as BT Men 

Smith, Dhillon fc Marsh 19981 , AC Cnc and V363 Aur 

Thoroughgood et al. 2004^ 



3.3 Light curves 

Fig. m shows the broad-band and emission-line light curves 
of V347 Pup. The emission-line hght curves were produced 
by subtracting a polynomial fit to the continuum and sum- 
ming the residual flux. All hght curves are plotted as a func- 
tion of phase following the ephemeris derived in Section r3.ll 
The b, y, V and i?-band hght curves show deep, 



asymmetrical primary eclipses with the egress lasting longer 
than ingress. Flickering is present in all hght curves, as 
well as an increase in brightness approaching eclipse in 
the b, y and l/-bands. The b and j/-band data recorded in 
1998 Jan show no significant brightness variations during 
the run, with out-of-eclipse magnitudes of 13.3 ± 0.1 
in both fllters. The eclipse depths are 3.2 and 2.6 mag, 
respectively. We measure _R-band out-of-eclipse magnitudes 
of 14.00 ± 0.05 mag on Dec 25, increasing in brightness to 
13.45 ± 0.10 mag on Dec 27. The eclipse depth remains 
roughly the same at 2.1 mag and 2.0 mag, respectively. 
In the l/-band, the out-of-eclipse magnitude is 14.10 ± 
0.10 mag, with an eclipse depth of 2.6 mag. Photometric 
out-of-eclipse magnitudes in the literature range between 
13.05-13.28 in R and 13.2-13.58 in V jBuckley et al. 1990| 
IMauche et al. 199411 . suggesting that our observations 
in Dec 1998 find V347 Pup around 0.5-1 mag fainter. 
Long-term variations in the magnitudes of NLs are not un- 
common (e.g. [Honeycutt 200l| l and have been observed in 
other SW Sex stars (e.g. BH Lyn, IDhillon et al. 19921 
DW UMa, IDhillon, Jones fc Marsh 1994| PX And, 
[Still, Dhillon fc Jone^^95f ! Low states are often ac- 
companied by the weakening or disappearance of the 
high-excitation Hell and CIII/NIII lines, which were 
unfortunately not observed in Dec 1998. There is, however, 
a change in the Hel A6678A Doppler maps between the 
two nights' observations, which is considered in Section f3.4l 
Further evidence that V347 Pup exhibits changes of 
state is seen in the EWs of the emission lines between 
the observed epochs. For example, the EW of H/3 varies 
between 17.0 ± 0.6A (July 1986, [Buckley et al. 1990l l, 62.6 
± 1.9 A (AprU 1991,lMauche et al. 1994II . 9.8 ± O.IA (Jan 
1995, [Still, Buckley fc Garlick 1998[ l and 24.6 ± 0.2A (Jan 
1998, this paper), although the high-excitation CIII/NIII 
complex has a constant EW between epochs. 

We measured the phase half-width of eclipse at the out- 
of-eclipse level (A0) by timing the first and last contacts of 
the eclipse and dividing by two. Our average value of Acj) 
= 0.110 ± 0.00 5 is consistent with the values of 0.120 ± 
0.011 quoted by jHarrop-Allin fc Warner (199 6)'and 0.105 ± 
0.005 measured by jBuckley et al. (1990)| We then computed 
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J-/ 1-L 



4200 4600 5000 5400 6200 6600 

Wavelength (A) 

Figure 1. The average spectra for the three wavelength regions; the spectrum centred on A4610A is an average of all spectra recorded 
on the Jan 1998 run, and has not been corrected for slit-losses. The spectrum centred on A5290A is an average of all data recorded on 
28 Dec 1998 and 23 Jan 1999, placed on an absolute flux scale (as determined from the 28 Dec 1998 photometry and flux standards). 
The spectrum centred on A6330A is composed of all data from 25 and 27 Dec 1998 and has been placed on an absolute flux scale. All 
average spectra are uncorrected for orbital motion, resulting in smeared spectral features. 



Table 3. Fluxes and widths of prominent lines in V347 Pup, measured from the 2 nights' data centred on A6330A and the night of 11 Jan 
1998 centred on A4610A. The full-width half-maximum (FWHM) velocities were determined from Gaussian fits, whereas the full-width 
zero-intensity (FWZI) velocities and their errors have been estimated by eye. Hell A4686A, CIII/NIII AA4640-4650A and Hel A4713A 
are blended, so separate values of the flux and EW are given (determined from a triple— Gaussian fit) as well as the combined flux of the 
three. 



Line Date Flux EW FWHM FWZI 

(x 10"" (A) (kms-i) (kms-i) 

(ergs cm~^ 



Ho 


25 Dec 


1998 


16.80 


± 


0.04 


35.1 


± 


0.1 


1100 


± 


100 


3500 


± 


300 


Ha 


27 Dec 


1998 


27.98 


± 


0.05 


36.4 


± 


0.1 


1100 


± 


100 


3600 


± 


300 


H/3 


11 Jan 


1998 


42.2 


± 


0.1 


24.6 


± 


0.2 


1000 


± 


100 


2800 


± 


300 


H7 


11 Jan 


1998 


34.6 


± 


0.2 


16.9 


± 


0.3 


1100 


± 


100 


2600 


± 


800 


Hel A4471A 


11 Jan 


1998 


8.0 


± 


0.1 


3.8 


± 


0.2 


1150 


± 


100 


1850 


± 


200 


Hel A4921A 


11 Jan 


1998 


4.23 


± 


0.08 


2.6 


± 


0.1 


1300 


± 


100 


2000 


± 


200 


Hel A5015A 


11 Jan 


1998 


3.3 


± 


0.1 


2.2 


± 


0.2 


1250 


± 


100 


2000 


± 


200 


Hel A6678A 


25 Dec 


1998 


1.49 


± 


0.02 


2.92 


± 


0.08 


1250 


± 


100 


1900 


± 


200 


Hel A6678A 


27 Dec 


1998 


2.40 


± 


0.04 


3.00 


± 


0.08 


1300 


± 


100 


1900 


± 


200 


CII A4267A 


11 Jan 


1998 


3.8 


± 


0.1 


1.7 


± 


0.4 


900 


± 


200 


1800 


± 


600 


Hell A4686A 


11 Jan 


1998 


26.9 


± 


0.3 


13.1 


± 


0.3 


1450 


± 


150 








CIII/NIII AA4640-4650A 


11 Jan 


1998 


11.7 


± 


0.1 


6.4 


± 


0.2 


1700 


± 


150 








Hel A4713A 


11 Jan 


1998 


4.8 


± 


0.3 


2.3 


± 


0.3 


1500 


± 


300 








Hell + CIII/NIII + 


11 Jan 


1998 


45.5 


± 


0.2 


22.4 


± 


0.2 















Hel A4713A 



the radius of the accretion disc in V347 Pup using the ge- 
ometric method outlined in |Dhillon, Marsh fc Jones (1991)] 
Combining A<j) with the system mass ratio and inclination 
derived in Section f3.1UI gives an accretion disc radius (Rd) 
of 0.72 ± 0.09 Ri, where Ri is the volume radius of the pri- 
mary's Roche lobe. This value is in agreement with the value 
of Rd/Ri > 0.82 quoted by |Harrop-Allin fc Warner (1996)| 
at the 2(7 level. 



The Ha eclipses are similar in shape to the continuum 
light curves, but do not appear to be as deeply eclipsed. 
The H/3 and H7 lines exhibit asymmetric eclipses, with 
ingress longer than egress. This behaviour is expected from 
asymmetric disc emission, consistent with the spiral arms 
identified in the Doppler maps (Section I8.4|l . The high- 
excitation Hell -I- CHI/NHI complex has a deep and U- 
shaped eclipse, suggesting an origin close to the white dwarf. 
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The Hel eclipses are wide with V-shaped minima, similar 
to the SW Sex stars (e.g. |Knigge et al. 2000l l. Note that the 
Hel flux is completely eclipsed, indicating an origin in the 
central portion of the disc, and not in an extended emis- 
sion region which is larger than the secondary star. The Hel 
A4471A emission line shows a broad dip in flux around phase 
0.4, before climbing to reach a maximum around phase 0.75, 
which could be a further signature of the disc asymmetry. 

3.4 Trailed spectrum & Doppler tomography 

We subtracted polynomial fits to the continuum and then 
rebinned the spectra onto a constant velocity-interval scale 
centred on the rest wavelength of the principal emission 
lines. For the data obtained in Jan 1998, we phase-binned 
all the spectra in order to boost the signal-to-noise (S/N). 
Individual spectra were weighted according to their S/N in 
order to optimally combine the spectra. The trailed spectra 
of Ha, H/?, Hell A4686A and Hel AA4471A, 5015A, 6678A 
are shown in Fig. |21 Doppler maps were calculated for the 
principal emission lines using the modulation Doppler to- 
mography code of [Steeghs (2003)| This method is an ex- 
tension to the conventional Doppler tomography technique 
(e.g. Marsh 2001 1, and maps both the constant and variable 
part of the line emission using a maximum-entropy regu- 
larised fitting procedure | |Skilling fc Bryan 1984} . We found 
that the modulated contribution to the line emission was 
weak (< 1 per cent), and thus our S/N was not sufficient 
to detect significant modulation in the accretion disc emis- 
sion. We therefore plot in Fig. 01 the corresponding aver- 
age Doppler maps only. The reconstructed line profiles are 
plotted next to the observed ones in Fig. |3 for comparison. 
Good fits to the data were achieved in all cases (reduced 
X' = 1 - 1-4) 

The Balmer-line trailed spectra are dominated by a 
low-velocity component with a semi-amplitude of ~150 
kms~^, moving from blue to red across primary eclipse. This 
is consistent with emission from the irradiated inner face of 
the secondary star, which is clearly seen in the correspond- 
ing Doppler maps. In the H/3 map, a second low-velocity 
emission source is present, seemingly coincident with the 
gas stream at a distance of 0.9Li, where Li is the distance 
from the white dwarf to the inner-Lagrangian point. There 
is also a weak two-armed disc asymmetry visible in the 
H/3 emission, which is much more prominent in the double- 
peaked Hel emission lines. Doppler maps of V347 Pup have 
been produced by [Still, Buckley fc Garlick (1998)| for data 
sets recorded in 1987, 1988 and 1995. The two compo- 
nents described above from the secondary star and the disc 
are clearly visible in their maps. The summed H/3 and H7 
maps of |Stm, Buckley fc Garhck (1998) [ show a stronger disc 
emission and spiral structure than our Balmer-line maps. 
The disc asymmetry is significant and is reminiscent of the 
two armed spiral structures that have been observed in the 
discs of dwarf novae during outburst (e.g. [Steeghs 2001| l. 
We return to these in Section [4. II The high-excitation Hell 
A4686A line is dominated by emission from the gas stream 
and bright spot overlayed on a weak accretion disc with 
radius Rd ~ 0.3 - 0.4Li. Note that the Hell A4686A 
Doppler map shows emission at higher velocities than the 
low-excitation lines, demonstrating that the material origi- 
nates from closer to the white dwarf. 



The blue and red Hel emission lines were recorded al- 
most a year apart and exhibit clear differences in structure. 
The secondary star emission is clearly evident in Hel A6678A 
(Dec 1998), afihough no strong Hel A447lA or A5015A emis- 
sion can be seen in the Jan 1998 data. There is also a dif- 
ference in the Hel A6678A Doppler maps between the 25th 
and 27th Dec 1998, which is probably related to the change 
in brightness of the system; the 25th Dec 1998 Doppler map 
has more enhanced spiral features and weaker secondary star 
emission than the 27th Dec (note that the average map of 
these two nights is shown in Fig.^J- During all these epochs, 
however, the spiral structures were observed, demonstrating 
that they are a persistent feature. 

3.5 Radial velocity of the white dwarf 

We measured the radial velocities of the emission lines 
in V347 Pup by applying the double-Gaussian method of 
Schneider fc Young (1980)] since this technique is sensitive 
mainly to the line wings and should therefore reflect the mo- 
tion of the white dwarf with the highest reliability. We tried 
Gaussians of widths 200, 300 and 400 kms~^ and we varied 
their separation from 200 to 3200 kms~^. We then fitted 

V = 7-iS'sin[27r(0-(/>o)] (2) 

to each set of measurements, where V is the radial ve- 
locity, K the semi-amplitude, (p the orbital phase, and 
(j)o is the phase at which the radial velocity curve crosses 
from red to blue. Examples of the radial velocity curves 
measured for the Ha, H^, Hell A4686A and Hel A447lA 
emission lines are shown in Fig. [Sj There is clear evidence 
of rotational disturbance in the emission lines, where 
the radial velocities measured just prior to eclipse are 
skewed to the red, and those measured after eclipse are 
skewed to the blue. This confirms the detection of a similar 
feature in the trailed spectra, and indicates that at least 
some of the emission must originate in the disc. There is 
also evidence of a phase shift in Ha and Hell A4686A, 
where the spectroscopic conjunction of each line occurs 
after photometric mid-eclipse. This phase shift implies 
an emission-line source trailing the accretion disc, such 
as a bright spot, and is a common feature of SW Sex 
stars (e.g. DW UMa, Shaffer, Hessman fc Zhan g 1988| 
V1315 Aql, [DhiUon, Marsh fc Jones 1991. SW Sex, 
[Dhillon, Marsh fc Jones 1997} . There appear to be no 
significant phase shifts, however, in the other emission 
lines. [Buckley et al. (1990)[ [Mauche et al. (1994)[ and 
Diaz fc Hubeny (1999)[ find no evidence of pheise shift in 
any of their emission lines, although their errors on (j)o were 
much larger. 

We tried to measure white dwarf radial ve- 
locity (Kw) values using a diagnostic diagram 
l |Shafter,~S zkody fc Thorstcnscn 1986), but with no 
success. We therefore attempted to make use of the 
light-centres method, as described by [Marsh (1988)1 the 
co-rotating co-ordinate system, the white dwarf has velocity 
(0, —Kw), and symmetric emission, say from a disc, would 
be centred at that point. By plotting Kx = —if sin (^0 
versus Ky = —Kcos(j)o for the different radial velocity 
fits (Fig. [SJ , one finds that the points move closer to the 
Ky axis with increasing Gaussian separation. A simple 
distortion which only affects low velocities, such as a bright 
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Figure 2. Broad-band and emission-line light curves of V347 Pup recorded in Dec 1998 and Jan 1999 (a), and Jan 1998 (b); see the 
panel labels for details. Note the increase in continuum and emission-line flux between Dec 25 and Dec 27 1998. 
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Figure 3. Trailed spectra and data computed from the Doppler maps (Fig. |4}. Tlie blue data recorded in Jan 1998 have been phase 
binned into 200s bins, the red data recorded in Dec 1998 into 300s bins. H7 has not been shown, as it is very similar in nature to H/3. 



spot, would result in this pattern, equivalent to a decrease 
in distortion as one measures emission further into the line 
wings and therefore more closely representing the velocity 
of the primary star. By linearly extrapolating the largest 
Gaussian separation on the Ha light-centre diagram to the 
Ky axis, we measure the radial velocity semi-amplitude of 
the white dwarf to be ~180 kms~^. The large uncertainty 
in this value (~ 40 kms~^), however, and the unsuccessful 



application of the technique to the other emission lines, 
prompted us to proceed with the mass determination using 
the secondary star features alone. 

3.6 Radial velocity of the secondary star 

The secondary star in V347 Pup is clearly visible in Fig. 
through absorption lines of Mgl, Fel and Cal. We compared 
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Figure 4. Doppler maps of the principal emission lines (H7 is not shown, as it is very similar in nature to H/3). The cross marked on 
each Doppler map represents the centre of mass of the system and the open circle represents the white dwarf. These symbols, the Roche 
lobe of the secondary star and the predicted trajectory of the gas stream, have been plotted using the /fij-corrected system parameters 
summarised in Table ISl The series of points along the gas stream mark the distance from the white dwarf at intervals of O.lLi, ranging 
from l.OLi at the red star to 0.2Li. Doppler tomography cannot properly account for variable line flux, so spectra around primary 
eclipse were omitted from the fits. 



regions of the spectra rich in absorption lines with a num- 
ber of templates with spectral types G7V-M3.5V. A tech- 
nique known as skew mapping was used to enhance the 
secondary features and obtain a measurement of the ra- 
dial velocity semi-amplitude of the secondary star (Kr). 
See |Vande Putte et al. (2003)|for a detailed critique of skew 
mapping and [Thoroughgood et al. (2004)] for a successful 
application to AC Cnc and V363 Aur. 

The data centred on A5290A were recorded specifically 
to exploit the secondary star features found between the H/3 
and Ha fines. Unfortunately, the presence of weak emission 
lines (e.g. Fell multiplet 42 at AA4924A, 5018A and 5169A, 
IMason fc Howell 20031 hampered all efi'orts to determine a 
Kr value from these data. The dominance of the emission 
lines in the spectra centred on A4610A also prevented a Kr 
determination from these data. The red spectra of V347 Pup 
centred on A6330A, however, allowed us to study the sec- 
ondary star through absorption features blueward of Ha, 
such as the Cal A6162A line. Exactly the same conclusion 
was reached by [Diaz fc Hubeny (1999)] 

The first step was to shift the spectral type template 
stars to correct for their radial velocities. We then normal- 



ized each spectrum by dividing by a constant and then sub- 
tracting a polynomial fit to the continuum. This ensures that 
line strength is preserved along the spectrum. The V347 Pup 
spectra were normalized in the same way. The template spec- 
tra were artificially broadened to account for both the or- 
bital smearing of the V347 Pup spectra due to their exposure 
times {texp), using the formula 

■ n2TTKR 



V = 



(3) 



(e.g. lWatson fc Dhillon 200111 . and the rotational velocity of 
the secondary (usini). Estimated values of Kr and vsini 
were used in the first instance, before iterating to find the 
best-fitting values given in Section r3.1UI Regions of the spec- 
trum devoid of emission lines were then cross-correlated 
with each of the templates, yielding a time series of cross- 
correlation functions (CCFs) for each template star. The 
regions used for the cross-correlation can be seen in Fig. |S] 
To produce the skew maps, these CCFs were then back- 
projected in the same way as time-resolved spectra in stan- 
dard Doppler tomography ijMarsh fc Horne 1988^ . If there 
is a detectable secondary star, we expect a peak at {0,Kr) 
in the skew map. This can be repeated for each of the tem- 
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Figure 5. Radial velocity curves of Ha, H/3, Hell A4686A and 
Hel A4471A using Gaussian widths of 300 kms~^ and a sepa- 
ration of 1400 kms~^. We omitted tlie points around primary 
eclipse during the fitting procedure (open circles) as these mea- 
surements are affected by the rotational disturbance. The emis- 
sion lines recorded in Jan 1998 have been phase-binned into 100 
bins for clarity. 




(km/s) 

Figure 6. Light-centres diagram for Ha. Points are plotted for 
radial velocity fits using Gaussians of FWHM = 300 kms-\ with 
separations from 900 kms~^ to 2900 kms~^ at 100 kms~^ inter- 
vals. The points move anti-clockwise, towards the Kx = axis 
with increasing Gaussian separation. 



plates, and the final skew map is the one that gives the 
strongest peak. 

The skew maps show well-defined peaks at Ky ~ 
216 kms~^ - the skew map of the M0.5V template is 
shown in Fig. |7] together with the trailed CCFs. A sys- 
temic velocity of 7 = 15 kms~^ was applied in order 
to shift the skew map peaks onto the Kx ~ axis 
(see [Smith, Dhillon fc Marsh 1998| for details). We there- 
fore adopt 7 = 15 ± 5 kms~ as the systemic veloc- 
ity of V347 Pup, in excellent agreement with the val- 
ues of 16 ± 10 kms~^ and 15 ± 12 kms^^ mea- 
sured bylStiU, Buckley fc Garlick (1998)] from the Balmer 
and Hell A4686l~ emission lines. The 7 velocities from 




-400 



-200 200 
Velocity (km/s) 



400 



Figure 7. Skew maps (top) and trailed CCFs (bottom) of 
V347 Pup cross-correlated with a M0.5V dwarf template. 



the emission lines shown in Fig. |3 ranged between 13 
kms^^ and 44 kms^^. Other 7 values measured from 
optical emission lines vary widely in the literature (- 
3 to 60 kms~\ [Diaz fc HubenyT999| -9 to 159 kms'^ 
IMauche et al. 199411 . 

Our adopted Kr of 216 ± 5 kms~^ was derived from 
the skew map peak of the best-fitting template found in Sec- 
tion E^Hl This result actually covers the Kh values derived 
from all of the template stars to within the errors, demon- 
strating that the result is robust to the choice of template 
(see Table lU. 



3.7 Rotational velocity and spectral type of the 
secondary star 

The spectral-type templates were broadened for smearing 
due to orbital motion, as before, and rotationally broadened 
by a range of velocities (50-240 kms~^). We then ran an op- 
timal subtraction routine, which subtracts a constant times 
the normalized template spectrum from the normalized aver- 
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Table 4. vsini values for V347 Pup cross-correlated with the 
rotationally broadened profiles of G7 - M3.5V templates. Also 
shown is the factor used to multiply the template star features 
during optimal extraction, and the position of the strongest peak 
in the skewmaps derived from each template using 7-velocities of 
kms^^ and 15 kms^^. 



Templates 


V sin i 


Optimal 








at min 


factor 


7 = 


7 = 15 




(kms-l) 




(kms"-"-) 


(km s^-"-) 


vjr / V 




n "^9 -t- n 






G9V 


133 


0.28 ±0.04 


(-13,215) 


(13,220) 


KOV 


133 


0.23 ± 0.03 


(-2,217) 


(14,219) 


KIV 


134 


0.24 ± 0.03 


(-15,215) 


(8,220) 


K2V 


133 


0.20 ±0.03 


(-22,212) 


(6,219) 


K3V 


136 


0.19 ±0.03 


(-28,212) 


(0,217) 


K4V 


135 


0.14 ±0.02 


(-17,211) 


(3,217) 


K5V 


134 


0.13 ±0.02 


(-17,213) 


(1,218) 


K7V 


133 


0.12 ±0.02 


(-24,210) 


(-3,216) 


M0.5V 


130 


0.13 ±0.02 


(-18,213) 


(0,216) 


M1.5V 


125 


0.12 ±0.02 


(-17,213) 


(-2,216) 


M2.5V 


126 


0.13 ±0.02 


(-21,213) 


(-7,216) 


M3.5V 


127 


0.12 ±0.02 


(-33,213) 


(-23,217) 



age V347 Pup spectrum, adjusting the constant to minimize 
the scatter in the residual. (Normalisation was carried out 
in the same way as Section 13.61 except that this time, the 
spectra were set to unity.) The scatter is measured by carry- 
ing out the subtraction and then computing the between 
the residual spectrum and a smoothed version of itself. By 
finding the value of rotational broadening that minimizes 
the 1 we obtain an estimate of both v sin i and the spec- 
tral type of the secondary star. Note that the v sin i values 
of the template stars are much lower than the instrumental 
resolution, so do not affect our measurements of wsini for 
the secondary star. 

The value of v sin i obtained using this method varies 
depending on the spectral type template, the wavelength 
region for optimal subtraction, the amount of smoothing 
of the residual spectrum in the calculation of the 
value of the limb-darkening coefficient used in the broaden- 
ing procedure. The values of vsini for all of the templates 
calculated using values for the limb-darkening coefficient of 
0.5 and smoothed using a Gaussian of FWHM = 15km s~^, 
are listed in Table 2] 

A plot of x^ versus v sin i for each spectral-type tem- 
plate is shown in Fig.lO] The spectral type with the lowest x^ 
value is M0.5V, which agrees with a visual identification of 
the best fitting template. [Diaz fc Hubeny (1999)| however, 
estimate a secondary star spectral type between KOV and 
K5V, with the possibility of a later-type subgiant. A plot of 
the V347 Pup average spectrum, a broadened M0.5V tem- 
plate spectrum and the residual of the optimal subtraction is 
shown in Fig.|H| The x'^ for the M0.5V template has a min- 
imum at 130 kms~^, so we adopt vsini = 130 ± 5 kms^^, 
with the error accounting for the measurement accuracy and 
the other variables noted in the previous paragraph. The er- 
ror quoted on our adopted value encompasses the measured 
V sin i for all of the templates used in the analysis (except 
for K3V with vsini = 136 kms~^). 



1.4 - 
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Figure 8. Orbitally— corrected average spectrum of V347 Pup 
(top) with the broadened M0.5V template (middle) and the resid- 
uals after optimal subtraction (bottom). The template spectrum 
has been multiplied by the scaling factor found from the opti- 
mal subtraction. All of the spectra are normalised and offset on 
the plot by an arbitrary amount for clarity. The wavelength lim- 
its shown are those used for the cross-correlation and optimal 
subtraction procedures, except for the region between the dashed 
lines owing to few secondary star features. 
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Figure 9. Determination of 11 sin i for V347 Pup using different 
spectral-type templates. Degrees of freedom = 699. 

3.8 The Kr correction 

The irradiation of the secondary stars in CVs by the 
emission regions around the white dwarf and the bright 
spot has been shown to influence the measured Kr (e.g. 
IWade fc Home 19881 I Watson fc Dhillon 200 l)l . For exam- 
ple, if absorption lines are quenched on the irradiated side 
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of the secondary, the centre of Ught will be shifted towards 
the back of the star. The measured Kr will then be larger 
than the true (dynamical) value. 

Diaz fc Hubeny (1999) [ found evidence for irradiation of 
the secondary star in V347 Pup, leading them to apply a 
correction to their measured Kr value. This fact, and the 
presence of Balmer and Hel emission from the inner face 
of the secondary star seen in the Doppler maps and trailed 
spectra fSection l3.4l l. prompted us to look for similar irradi- 
ation effects in the absorption lines of our data. We applied 
the following two observational tests. First, the rotation- 
ally broadened line profile would be distorted if there was 
a non-uniform absorption distribution across the surface of 
the secondary star ( |Davey fc Smith 1992[ l. This would re- 
sult in a non-sinusoidal radial velocity curve. Second, one 
would expect a depletion of secondary star absorption-line 
flux at phase 0.5, where the quenched inner-hemisphere is 
pointed towards the observer (e.g. lFriend et al. 1990t . 

The secondary star radial velocity curves were pro- 
duced by cross-correlating the V347 Pup spectra with the 
best-fitting smeared and broadened template spectra, as 
described in Section 13.61 The cross-correlation peaks were 
plotted against phase to produce the radial velocity curves 
shown in the lower panel of Fig. 1101 There is evidence for an 
eccentricity in the radial velocity curve compared with the 
sinusoidal fit represented by the thin solid line, although the 
data are noisy. 

The variation of secondary star absorption-line flux 
with phase for V347 Pup is shown in the top panel of 
Fig. 1101 These light curves were produced by optimally 
subtracting the smeared and rotationally broadened best- 
fitting template from the individual CV spectra (with the 
secondary radial velocity shifted out) as described in Sec- 
tion 13.71 This time, however, the spectra were continuum 
subtracted rather than normalised to ensure that the mea- 
surements were not affected by a fluctuating disc bright- 
ness. The constants produced by the optimal subtraction 
are secondary star absorption-line fluxes, correct relative to 
each other, but not in an absolute sense. The dashed lines 
super-imposed on the light curves represent the variation 
of flux with phase for a Roche lobe with a uniform absorp- 
tion distribution. The sinusoidal nature is the result of the 
changing projected area of the Roche lobe through the or- 
bit. The V347 Pup light curve is clearly not represented 
by a uniform Roche lobe distribution as the secondary star 
absorption-line flux vanishes between phases 0.4-0.6. 

These three pieces of evidence, as well as the disap- 
pearance of the CCFs between phases 0.4-0.6 seen in Fig.[7| 
suggest that the secondary star in V347 Pup is irradiated 
and we must correct the Kr values accordingly. 

It is possible to correct Kr for the effects of irradia- 
tion by modelling the secondary star flux distribution. In 
our simple model, we divided the secondary Roche lobe 
into 40 vertical slices of equal width from the Li point to 
the back of the star. We then produced a series of model 
light curves (using the system parameters derived in Sec- 
tion 13.101 , varying the numbers of slices omitted from the 
inner hemisphere of the secondary which contribute to the 
total flux. The model light curves were then scaled to match 
the observed data, and the best-fltting model found by mea- 
suring the between the two. In all models, we used a 
gravity-darkening parameter /3 = 0.08 and limb-darkening 



coeflicient u = 0.5 fe.g. IWatson fc Dhillon 200H . The neg- 
ative data points around phase 0.5 were set to zero, as the 
secondary star absorption line flux disappears at this point. 
Once the best-fitting light curve was found, we produced 
fake V347 Pup spectra from the model, which were cross- 
correlated with a fake template star to produce a synthetic 
radial velocity curve. In the first instance, the synthetic 
curve mimicked the non-sinusoidal nature of the observed 
data, but with a larger semi-amplitude. This was expected, 
as the model input parameters used the uncorrected Kr de- 
rived in Section r3.10l We then lowered Kr and repeated the 
process, until the semi-amplitude of the model and observed 
radial velocity curves were in agreement, each time checking 
the light curve models for goodness of fit. The resulting Kr 
was then adopted as the real (or dynamical) Kr value. 

The best-fltting model light curve was produced by 
omitting 12 slices when fitting the data (reduced between 
model and data = 1.03). The model light curves omitting 11, 
12 and 13 slices are shown by the solid lines in Fig. 1101 Our 
final model, which has an input Kr of 198 kms~^, produces 
the radial velocity curve shown as the thick solid line in the 
lower panel of Fig. 1101 There is good agreement between 
this and the observed data. If gravity-darkening and limb- 
darkening are neglected, the best flt light curve remains the 
same, but produce a Kr value which is ~ 6 kms~^ lower. 

In summary, we correct the Kr of V347 Pup from 216 
kms~^ to 198 kms~^. This correction of 18 kms"'^ is ex- 
actly the same as that calculated by |Diaz fc Hubeny (1999)] 
using a much simpler approximation, which changed their 
measured value of 205 kms~^ to 187 kms~^. 



3.9 The distance to V347 Pup 

By finding the apparent magnitude of the secondary star 
from its contribution to the total light of the system, and 
estimating its absolute magnitude, we can calculate the dis- 
tance (d) using the equation: 

5 \og{d / 10) ^mv - Mv - dAv/ WOO (4) 

where Av is the visual interstellar extinction in magnitudes 
per kpc. 

The mean i?-band photometric flux of V347 Pup during 
the recorded spectra is 8.93 mjy, which we convert to a mean 
7?-band magnitude of 13.8 ± 0.3. The uncertainty reflects 
the change in brightness of the system between 25 and 27 
Dec. During this time, the secondary star contributes 13 ± 
2 per cent of the total light of the system, assuming an early 
M spectral type (see Table QJ. The apparent magnitude of 
the secondary is therefore R = 16.0 ± 0.4, which we convert 
to a V^-band magnitude of 17.3 ± 0.4 using a typical V — R 
value for an early M star from Gray (1992)] There are a 
number of ways of estimating the absolute magnitude of the 
secondary star, assuming it is on the main sequence (e.g. 
IPatterson 19841 IWarner 199*5^1 [Gray 1992t . We took each 
of these into account and adopt an average value of Mv ~ 
+8.8 ± 0.5. [Mauche et al. (1994)| estimated the extinction 
to V 347 Pup to be E {B -V) = 0.05, which results in Ay = 
0.16 HScheffl er 198211. The distance to V347 Pup is calculated 
from equation 01 to be 490 ± 130 pc. 

Buckley et al. (1990)| estimate the distance to 
V347 Pup to be between 174-380 pc, based on their 
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Orbital phase 

Figure 10. Upper panel: Secondary star absorption line light 
curve with model fits (solid lines). Model fits are shown for Roche 
lobes with 11, 12 and 13 slices removed (see text for details). The 
lower the line, the more slices removed. The dashed line represents 
a model where slices have been removed. The data have been 
phase— binned into 50 bins to increase S/N. Lower panel: Measured 
secondary star radial velocity curve with a sinusoidal fit (thin 
solid line) and the best— fitting model fit (thick solid line). In both 
panels, the open circles indicate points that have been omitted 
from the fits and the data have been folded to show 2 orbital 
phases. 



measured system inclination and out-of-eclipse magni- 
tude. [Mauche et al. (1994)| use their interstellar reddening 
measurement and a mean interstellar hydrogen number 
density to estimate a distance of 340-590 pc. Finally, 
[Diaz fc Hubeny (1999)| find a distance of 510 ± 160 pc from 
the spectral line depths of the secondary star. Our value is 
consistent with all distance estimates in the literature. 



3.10 System parameters 

Using the Kr and v sin i values found in Sections l8.7l and l8.8l 
in conjunction with the period determined in Section 13.11 
and a measurement of the eclipse full-width at half depth 
(A01/2), we can calculate accurate system parameters for 
V347 Pup. 

In order to determine l\4>i/2, we estimated the flux out 
of eclipse (the principal source of error) and at eclipse mini- 
mum, and then measured the full-width of the eclipse half- 
way between these points. The eclipse full-width at half- 



depth was measured to be A0i/2 = 0.115 ± 0.005, in agree- 
ment with the eclipse half-width at half-depth of 0.052 ± 
0.002 measured by [Buckley et al. (1990)] at the 2a level. 

We have opted for a Monte Carlo approach similar to 
[Home, Welsh fc Wade (1993)| to calculate the system pa- 
rameters and their errors. For a given set of Kr, vsini, 
/S.(f>i/2 and P, the other system parameters are calculated 
as follows. 

i?2/a can be estimated because we know that the sec- 
ondary star fills its Roche lobe (as there is an accretion disc 
present and hence mass transfer) . R2 is the equatorial radius 
of the secondary star and a is the binary separation. We used 
Eggleton's formula | |Eggleto~ 1983J which gives the volume- 
equivalent radius of the Roche lobe to better than 1 per cent, 
which is close to the equatorial radius of the secondary star 
as seen during eclipse. 



R2 



0.49g- 



2/3 



(5) 



a 0.6g2/3 +ln(l + gi/3)- 

The secondary star rotates synchronously with the orbital 
motion, so we can combine Kr and vsmi, to get 

— (i + a) = -^- (6) 

a Kr 

By considering the geometry of a point eclipse by a spherical 
body (e.g. Dhillon, Marsh fc Jones 1991| l, the radius of the 
secondary can be shown to be 



] — siv? 7rA(;/>i/2 -f cos^ 7rA0i/2 cos^ i, 
a ' 



(7) 



which, using the value of R2/a obtained using equations |3 
and|S| allows us to calculate the inclination, i, of the system. 
The geometry of a disc eclipse can be approximated to a 
point eclipse if the light distribution around the white dwarf 
is axi-symmetric (e.g. IDhillon 1990t . This approximation is 
justified given the symmetry of the primary eclipses in the 
photometry light curves (Fig. |2J. Kepler's Third Law gives 
us 



KlPorb _ Ml sin'^ 

2^G " (l + g)2 



(8) 



which, with the values of q and i calculated using equations 
El and 13 gives the mass of the primary star. The mass of 
the secondary star can then be obtained using 



M2 = qMi 



(9) 



The radius of the secondary star is obtained from the equa- 
tion 

V sin i 2n sin i 

(e.g. IWarner 1995all and the separation of the components, 
a, is calculated from equations |S|and llOl with q and i now 
known. 

The Monte Carlo simulation takes 10 000 values of Kr, 
vsini, and A(^i/2 (the error on the period is deemed to be 
negligible in comparison to the errors on Kr, vsini, and 
^01/2), treating each as being normally distributed about 
their measured values with standard deviations equal to the 
errors on the measurements. We then calculate the masses of 
the components, the inclination of the system, the radius of 
the secondary star, and the separation of the components, as 
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Figure 11. Monte Carlo determination of system parameters for 
V347 Pup. Each dot represents an Mi , M2 pair; the solid curves 
satisfy the vsini and Kji constraints, and the dashed lines mark 
lines of constant inclinations (i = 80°, 85° and 90°). 



outlined above, omitting (ifa, usini, ^(^y^j^) triplets which 
are inconsistent with sini < 1. Each accepted Mi,M2 pair 
is then plotted as a point in Figure fTTI and the masses and 
their errors are computed from the mean and standard de- 
viation of the distribution of these pairs. 

We find the component masses of V347 Pup to be 
Ml = 0.63 ± O.O4M0 and M2 = 0.52 ± 0.06Mq. The values 
of all the system parameters deduced from the Monte Carlo 
computation are listed in Table |H1 including i^^ij-corrected 
and non /S'^-corrected values for comparison. Note that our 
derived Kw of 163 ± 9 kms~^ is in remarkable agreement 
with the Kw values of |Still, Buckley fc Garlick (1998)1 who 
measure 156 ± 10 kms~ using a double-Gaussian convo- 
lution of the Balmer lines, and 166 kms~^ as the centre 
of axisymmetric Balmer emission. The white dwarf mass 
of 0.63 ± O.O4M0 is consistent with the average value of 
Ml = 0.80 ± O.22M0 (for CVs abov e the period gap) de- 
termined bylSmith & Dhillon (1998)1 The empirical relation 
obtained by Smith fc Dhillon ( 1998) | between mass and ra- 
dius for the secondary stars in CVs predicts that if the sec- 
ondary star in V347 Pup is on the main-sequence, it should 
have a radius of 0.54 ± 0.08 Rq. Our measured value of 0.60 
± 0.02 Rq (from equation 1101 is consistent with this value. 



4 DISCUSSION 
4.1 Spiral arms 

Spiral-armed disc asymmetries are evident in the 
Hel and H/3 Doppler maps, confirming the find- 
ings of [Still, Buckley fc Garlick (1998)| in their H/3 
and H7 maps. Similar spiral structures have been 
observed in dwarf novae during outburst (e.g. 
IP Peg, [Steeghs, H arlattis fc Ho me 1997| U Gem, 
[Groot, Rutten i^TTarTT'aradijs 2001 J . Tidally driven 



spiral density waves can develop in accretion discs due 
to the tidal torque of the mass donor star on the outer 
disc (ISawada, Matsuda fc Hachisu 1986| IBTondin 20Qni 
IBoftin 20^r Their detection in outburst only reflects 
the much stronger tidal effects on the accretion disc 
when it increases in size and temperature during out- 
burst, in which case a tidally induced spiral structure 
is expected that closely matches the observed struc- 
tures (|Armitage fc Murray 1998| [Steeghs fc Stehle 1999| 
[Steeghs 2001| |. in d-irarf novae, these asymmetries decay 
as the system returns to quiescence, and the disc cools 
and shrinks. In order for a similar tidal response to be 
responsible for the disc asymmetry in V347 Pup, its disc 
must be large and comparable to the tidal radius. We 
calculate the tidal radius of the accretion disc to be 0.33a 
using the pressureless disc models of [Paczyriski (1977) [ and 
our new system parameters. The measured disc radius of 
Rd/cl = 0.28 ± 0.03 is comparable in size to the tidal 
radius, and therefore consistent with a tidal origin for the 
observed spiral structure. 

Our observations show that the spiral structures are 
clearly visible in the Hel Doppler maps, but are either weak 
or non-existent in the Balmer and Hell maps. This is in 
contrast to dwarf novae in outburst, which typically show 
stronger spiral structures in the Hell and Balmer lines (e.g. 
[Mars h 2 001 Morales- Rueda 2_004j. This could be a reflec- 
tion of different densities and temperatures in NL discs com- 
pared to the discs of dwarf novae in outburst, or it could 
simply be due to a contrast effect where the relative contri- 
bution of the spiral structure is not as high in the Hell and 
Balmer maps due to the presence of low-velocity emission. 

Note that the impact of such tidally-induced spi- 
ral arms on the angular momentum transport has not 
been fully established. If they are associated with hy- 
drodynamical shocks, such as in the simulations of 
[Sawada, Matsuda fc Hachisu (1986)[ their contribution to 
the angular momentum transport could be very significant. 
On the other hand, [Smak (2001)[ and [Ogilvie (2002) [ propose 
that these disc structures may reflect tidally thickened ar- 
eas in the outer disc as it expands close to its tidal radius. 
Their enhanced emission is then caused by irradiation from 
the accreting white dwarf and regions close to it. 

The prospect of testing such basic disc physics with 
observations warrants the study of these disc structures in 
more detaU (see also lMorales-Rueda 200411 . With V347 Pup, 
we have a target that appears to have a persistent disc asym- 
metry that is more accessible than the transient spiral struc- 
ture observed in dwarf novae. 



4.2 Mass transfer stability 

The mass ratio of a CV is of great significance, as it governs 
the properties of mass transfer from the secondary to the 
white dwarf primary. This in turn governs the evolution and 
behaviour of the system. 

The secondary star responds to mass loss on two 
timescales. First, the star returns to hydrostatic equilibrium 
on the dynamical timescale, which is the sound-crossing 
time of the region affected. Second, the star settles into a new 
thermal equilibrium configuration on a thermal timescale. 

The two timescales upon which the secondary responds 
to mass loss leads to two types of mass transfer instability. 
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Table 5. System parameters for V347 Pup. The Monte Carlo results for corrected and uncorrected Kji values are shown for comparison. 
The radial velocity of the white dwarf {Kw) has also been calculated from the secondary star parameters. 
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If, upon mass loss, the dynamical response of the secondary 
is to expand relative to the Roche lobe, mass transfer is 
dynamically unstable and mass transfer proceeds on the dy- 
namical times cale. [Politano (1996)] made an analytic fit to 
the models of |Hjellming (1989)| to give the limit of dynam- 
ically stable mass loss, plotted as the solid line in Fig. 1121 
Dynamically stable mass transfer can occur if the CV lies be- 
low this line. This limit is important for low mass secondary 
stars (M2 < 0.5Mq), as they have significant convective 
envelopes that tend to expand adiabatically in response to 
mass loss 1^ Kool 19921. 

Thermally unstable mass transfer is possible if the dy- 
namic response of the star to mass loss is to shrink rela- 
tive to its Roche lobe (i.e. mass transfer is dynamically sta- 
ble). This occurs at high donor masses (M2 > O.8M0) when 
the star has a negligible convective envelope and its adia- 
batic response to mass loss is to shrink, (e.g. Ide Kool 19921 
IPolitano 1996t . Mass transfer then initially breaks contact 
and the star begins to settle into its new thermal equilibrium 
configuration. If the stars thermal equilibrium radius is now 
bigger than the Roche lobe, mass transfer is again unstable, 
but proceeds on the slower, thermal timescale. The limit of 
thermally-stable mass transfer can be found by differenti- 
ating the main-sequence mass-radius relationship given in 
[PolitanoJJl996)| Thermally-stable mass transfer can occur 
if the CV appears below the dotted line plotted in Fig. 1121 

The limit for dynamically stable mass transfer is im- 
portant in the case of V347 Pup owing to the low secondary 
star mass. Fig. 1121 shows that the system is just consistent 
with the limit at the la level. The mass transfer stabil- 
ity limits, however, are true only for ZAMS stars, whereas 
the secondary stars in CVs are expected to have undergone 
some evolution. The loss of the outer envelope, for example, 
would result in a larger than normal helium to hydrogen ra- 
tio and affect the star's response to mass loss. For instance, 
DX And, which lies outside the limit, has been shown to 
have an evolved companion ( |Drew, Jones fc Woods 19931 - 



There is tentative evidence that the secondary star 
in V347 Pup is evolved by considering three pieces of ev- 
idence. First, V347 Pup falls outside the limit for dy- 
namically stable mass transfer (although agrees at the 
la level). Second, the measured radius is at the upper 
limit for a main-sequence companion of the same mass 
llSmith fc Dhillon 199 81. Third, the secondary star mass 
and spectral type measured for V347 Pup are closer to 
the evolved models of |Kolb, King fc Baraffe (2001)| than the 
ZAMS models. 



5 CONCLUSIONS 

(i) We have measured the radial and rotational velocities 
of the secondary star in V347 Pup in order to calculate the 
component masses and other system parameters. The sec- 
ondary star radial velocity is affected by irradiation from the 
emission regions around the primary, which we correct for 
using a model. We find the component masses in V347 Pup 
to be Ml = 0.63 ± 0.04 Mq for the white dwarf primary 
and M2 = 0.52 ± 0.06 Mq for the M0.5V secondary star. 
V347 Pup shows many of the characteristics of the SW Sex 
stars, exhibiting single-peaked emission lines, high-velocity 
S-wave components and phase-offsets in the radial velocity 
curves. 

(ii) V347 Pup lies outside the theoretical limit for dynam- 
ically stable mass transfer in ZAMS stars, but is just con- 
sistent at the la uncertainty level. This piece of evidence, 
together with a secondary star radius at the upper limit 
for a main-sequence companion of the same mass, suggests 
that the secondary star in V347 Pup may be evolved. Addi- 
tionally, the secondary star mass and spectral type mea- 
sured for V347 Pup are closer to the evolved models of 
|Kolb, King fc Baraffe (200T]1 than the ZAMS models. 

(iii) The presence of spiral arms in the accretion disc, 
first noted by |Still, Buckley fc Garlick (1998)| has been con- 
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Figure 12. Critical mass ratios for mass transfer stability. 
The dotted line represents the condition for thermal instabil- 
ity; the solid line represents the condition for dynamical insta- 
bility (Politano 1996). Both curves assume the star is initially 
in thermal equilibrium. Mass ratios and secondary masses from 
the c ompilation of I Smith fc D hillon (1998)' 'North et al . (2000)] 
[Wats on ct al. (2003)| and [Thoroughgood ct al. (2004) are over- 
plotted. The mass ratios and secondary star masses of V347 Pup 
determined in this paper are also plotted. 



firmed. Consistent witli this, we find that the measured ac- 
cretion disc radius is close to the tidal radius computed from 
the pressureless disc models of jPaczynski (1977)| The per- 
sistent spiral arms seen in this bright novalike makes it an 
excellent candidate in which to study these features, rather 
than the transient spiral structures observed in dwarf novae. 
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